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ABSTRACT 



This report describes the procedure utilized in the assembly and 
evaluation of a small intermittent blowdown supersonic wind tunnel. 

The tunnel proper was fully constructed and nozzle blocks for 
attainment of a Mach number of 1.5 or 2.37 in the 1/2*’ by 1.718” 
test section were available. The basic schlieren system was purchased 
complete. An attempt was made to use or convert readily available 
materials for all components. 

It was expected the results of this project would be of a perman- 
ent nature and would be utilized by the Mechanical Engineering Depart- 

■* 

ment as a laboratory installation. The installation is suitable for 
the study of compressible fluid flow and the phenonemon of normal 
shocks but is not suitable for the study of the performance of shapes 
or models in the fluid flow. 

The tunnel installation operates satisfactorily with a maximum 
Mach number of 2.29 with approximately 0.28 Ibm of air per second. 
Schlieren photographs of the flow field may be taken in either black 
and white or in color. 
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INTRODUCTION 



There are two basic types of wind tunnels available for the study 
of compressible fluid flow (1) the continuous and (2) the intermittent 



pressor, tunnel volume flow and pressure ratio characteristics present 
major difficulties in the design work. In addition the Mach number span 
of a continuous wind tunnel is usually small and does not extend from 
low subsonic velocities to high supersonic velocities. These factors 
combined with a higher cost of installation and operation of the con- 
tinuous tunnel have resulted in a shift towards the use of the inter- 
mittent tunnels. 

The intermittent tunnel is further divided into two basic types^^^(l) 
the vacuum storage drive and (2) the pressure storage drive, [[^ll]] . In 
the first system, atmospheric air passes through a drier, wind tunnel 
and diffuses into a vacuum vessel which is continuously evacuated to the 
atmosphere. In the second system, dry air is stored under pressure in 
a storage tank; the air passes through a throttle valve, wind tunnel, 
and exhausts through a tunnel diffuser to the atmosphere. 

It was the objective of this thesis to install and evaluate the per- 
formance of a small pressure storage type supersonic wind tunnel system. 
The components available at the start of the project consisted oft high 
pressure compressor; storage tanks; two regulator valves; tunnel test 
section; and Schlieren system components; all other pieces for the 
system had to be either specified or designed and specially constructed. 
It was further necessary to carry out assembly and installation of these 
components. Having completed this phase of the objective, it was then 

* Numbers in brackets refer to the Bibliography on page 20. 




With continuous wind tunnels the problem of matching com- 
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possible to undertake a performance evaluation program. This con- 
sisted of both a qualitative and an quantative approach. The qualita- 
tive evaluation was derived primarily from Schlieren photographs^ while 
the quantitative results came primarily from pressure surveys. 

Figure 1 is a schematic diagram of the system showing the basic 
design and necessary components. Discussion of the design and per- 
tinent details of certain sections of the system is found in Part I. 

Figure 2 is a subassembly drawing of the small tunnel test section 
for which the performance tests were made. Results of these tests are 
presented in photographs, graphical form, and tabular form. Discussion 
of these tests and results appear in Part II. 

The Appendicies include various sample calculations, a discussion 
of Schlieren systems, operational instructions, design of a differen- 
tial pressure transducer, and design of needed heat capacitors for the 
air receivers. 
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Part I. Tunnel System Design 

A. Air Supply , 

Air for the system was supplied by a Chicago Pneumatic two stage 
inner cooled air compressor rated at 114 cubic feet of free air per min-^ 
ute or 0.142 pounds per second. The compressor is driven by a 50 horse- 
power induction electric motor. The compressor discharge pressure can 

be controlled up to 500 psi by the regulator valve (Bailey Type-XSDS-9M) . 

% 

Operational instructions for the compressor are given in Appendix A. 

Air from the compressor goes to an after-cooler (Graham Heliflow. 
Type 9AA 145) where the temperature is reduced to about 60°Fs and then 
to a cyclone type separator with trap (Wright Austin)^ where the con- 
densed water is removed. Calculations show the moisure content after 
the air leaves the separator to be about 0.000319 lb. water per lb of 
dry air, Appendix B. Drying of the air is necessary to prevent condensa- 
tion from occuring in the tunnel test section, \j>^ , Ij-ll s> • 

From the separator, the air passes through an oil filter (C. M. 

Kemp, Type FPF-1610), using activated alumina (Vapoilsorb) as a desi- 
ccant. This filter was modified by putting a 60 mesh stainless steel 
strainer on the inlet side along with a fiberpacked strainer. This modi- 
fication was made to prevent blowback of desiccant particles to the 
compressor by high pressure air when the compressor was unloaded. 

The air then goes through the regulator valve, which maintains 
a constant compressor discharge pressure, to the two air receivers. 

The receivers, connected in parallel, provide a total storage volume of 
234 cubic feet. Figure 3 is a photograph showing these various compon* 
ents in relation to the acoustic enclosure which houses the compressor. 



3 



When operating the wind tunnel, the mass flow of air from the 
receivers to the tunnel section is controlled by a regulator valve 
(Bailey, Type XSDS-9M). The regulator valve provides for stagnation 
pressure control up to 200 psi. Figure 4 shows the position and loca- 
tion of this valve. 

B. Tunnel Section . 

The tunnel section consists of six basic parts: stilling tank; 

entrance nozzle; test section; diffuser; flow meter; and back pressure 
control. 

The stilling tank was made from a 30 inch piece of nominal 6 inch 
pipe flanged at both ends. The flanges were faced and grooved for use 
with 0-ring seals. The upstream end of the tank was connected by 
means of a piece of high pressure rubber hose to the regulator valve 
discharge. Figure 5 is a photograph of these pieces. In the upstream 
end of the stilling tank, 3 stainless steel 60 mesh screens were used 
to diffuse any high velocity profile. 

The entrance nozzle was constructed of a polyester resin and 
chopped fiber glass molded over a wooden core, A pair of cubic equa- 
tions were used to define the nozzle contour so that flow separation 
would not occur. The flanges were molded integral with the nozzle and 
were faced and drilled after casting. The upstream flange has a 1/2 
inch lip for a matching connection with the stilling tank. 

The test section was designed and constructed by the Ames 
Aeronautical Laboratory, Moffett Field, California. Features of this 
test section consist of a variable area supersonic diffuser, inter- 
changeable nozzle blocks, and interchangeable side walls for either 
visual or pressure surveys. 
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The present sets of nozzle blocks are for Mach numbers 1.50 and 
2.37. The sidewall for pressure surveys is provided with 32 pressure 
taps. Figure 2 is a subassembly drawing of the section and Figure 6 
is a photograph showing the section with side walls removed. Tables I 
and II give the nozzle coordinates, area ratios, and design Mach numbers 
and pressure ratios based on one-dimensional calculations for the two 
sets of nozzle blocks. 

The diffuser, like the entrance nozzle, was also molded over a 
wooden core with a polyester resin and fiber glass. To provide added 
strength, 1/8 by 3/8 inch steel bars were molded into the side and 
flanges. The design was straight wall with a maximum divergance angle 
of 7® to prevent flow separation; and diffusion of the air stream was 
carried to a point where the flow could be treated as incompressible for 
metering purposes. 

The four pieces just described were integrally mounted to a movable 
stand and are shown in Figure 7. 

The flow meter was designed according to A.S.M. E. standards for a 
sharp edge orifice with D and (1/2)D pressure taps, W . fs] . An en- 
trance length of 20 pipe diameters and an exit length of 10 pipe diameters 
were provided. Nominal 3 inch pipe was used which was bored to a 3.120 
inch diameter for 16 inches either side of the orifice. 0-ring seals 
were used throughout including the flanged joint with the diffuser. 

Orifice plates were made of 1/8 inch stainless steel. Appendix C presents 
the details for computing mass flow rates. 

The back pressure control consisted of a manually operated 3 inch 
globe valve. This control was incorporated so that changes could be made 
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In the test section air density. Since this throttling was somewhat 
noisy, a large truck muffler was added to the valve exit® 

These pieces were also mounted on a movable stand and are shown 
in Figure 8. 

The completely assembled tunnel section is shown in the photo- 
graph, Figure 9. 

C. Instrumentation . 

Instrumentation for the tunnel section consisted of devices for 
three types of measurement; pressure, temperature, and velocity. 

Manometers were used for all quantitative pressure measurements 
while bourdon tube gages were lised for qualitative check purposes. 

Two 120 inch mercury differential manometers were utilized with the 
32 tap pressure survey sidewall of the test section. Two 16 connector 
manifolds with quick opening valves were made to accommodate the 32 taps 
to the two manometers. Figure 10 shows the arrangement of the manifolds 
with the manometers. It was originally planned to produce strain gage 
type differential pressure transducers for these measurements to de- 
crease the reading time required; fabrication difficulties, however, 
prevented the completion of this project, Appendix D. 

A 120 inch mercury differential manometer was also used to monitor 
the test section throat pressure continuously during pressure survey runs. 

Stagnation pressure, from the stilling tank tap, was likewise con- 
tinuously monitored with a 120 inch mercury manometer besides being fed 
to the differential manometers. 

Static pressure for the flow meter could be read with either an 
oil (sp. gr. = 0.834) or mercury filled 60 inch manometer. Differen- 
tial pressure for the orifice was determined with a water manometer. 
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Ellison inclined manometers were used in conjunction with the pitot 

tube. 

The only quantitative temperature measured was the stagnation 
temperature in the stilling tank. This was accomplished with a 
stainless steel radiation shielded 30 gage copper-constantan thermocouple 
referenced to ice. The output emf of the thermocouple was read with a 
continuous balance precision potentiometer (Minncapolis-Honeywellj» Model 
SY156X-VH). Various other temperatures of the tunnel system were meas- 
ured with gas bulb indicating thermometers. 

Velocity measurements in the stilling tank were made with a 1/8 
inch diameter Prandt 1-Pitot p^robe. These measurements were made sole- 
ly for the purpose of determining the uniformity of the velocity pro- 
file leaving the straightening screens in the stilling tank and were 
only of a qualitative nature. The Pitot probe was removed from the 
stilling tank after finding that the screens adequately diffused the 
high velocity entering air to a uniform low velocity profile. 

Velocity measurements and visual flow studies in the test sec- 
tion were accomplished with a coincidence type Schlieren systeii manu- 
factured by Tinsley Laboratories Inc. A schematic diagram is shown in 
Figure 11. The components were mounted to a rigid moveable stand with 
the wheels set in tracks parallel to the test section, Figure 8. The 
viewing section of the Schlieren system was modified by installation 
of an Ilex Synchro shutter and a 4 by 5 inch ground glass view plate 
and film pack holder. Figure 12. This modification allowed photographs 
of the flow behavior to be made which, when enlarged in printing^ pro- 
vided for qualitative and some quantitative measurements. Appendix E 
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presents a general theory of the Schlieren method of flow visualization 



as well as 



operating instructions for this particular system. 
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Part II, Performance Evaluation 



A. Qualitative Results , 

The qualitative evaluation of the tunnel system performance was 
derived primarily from Schlieren photographs. Figure 13 depicts the 
tunnel section operation under shock free conditions. The knife edge 
intercepted the light beam in both a vertical and horizontal position 
to show the change of density patterns in the horizontal and vertical 
directions. Figure 13c, with the knife edge horizontal, clearly shows 
the boundary layer. 

Criss-crossing the tunnel outline are Mach wavelets produced by 
wall friction and other irregularities of the tunnel surfaces. These 



wavelets are related 






to the Mach number of the air flow in the 



tunnel by the relationship, 
sin (?C - PI 

where <?C is the acute angle between the wavelets and a horizontal line. 
Figure 18a is an enlargement of a section of Figure 13a and Figure 18b 
is an enlargement of Figure 13b, The angle in Figure 18a is 27.50° 
with a resulting Mach number of 2.17. This wavelet originated at the 
junction of the nozzle and the test section. Utilizing the pressure 
ratios, the Mach number at the corresponding point was determined to 
be 2.125, The angle in Figure 18b is 29.25° with a resulting Mach 
number of 2,04. The Mach number calculated for this position was 1,828. 
This wavelet originated at approximately pressure tap 30. 

Figure 14 depicts a normal shock occuring in the diffuser section 
which was being held at a constant height of 1.718 inches. Inner action 
between the shock and the boundary layer is clearly indicated. Up- 
stream of the shock the boundary layer is stable. At station 12 there 
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is a separation of the boundary layer due to the interaction of the 
shock and the boundary layer. Less severe interaction occurs between 
the boundary layer and the repeated shocks do\^mstream until complete 
subsonic flow is established downstream of station 15. 

Comparing Figure 14a operating under a stagnation pressure of 
67 inches of mercury and a back pressure of 30.5 inches of mercury to 
Figure 14c for tunnel operation with a stagnation pressure of 81.5 inches 
of mercury and a back pressure of 38.3 inches of mercury, it may be seen 
that the resultant shock has moved upstream almost three inches. The 
pressure ratio at station 11 in Figure 14a v;here no shock exists is 
approximately 0.1088 as taken from Figure 23 or Table V. At a point 
forward of the shock in Figure 14c a pressure ratio of 0.103 exists 
as taken from Table VII. The increased density of the air and the high- 
er Reynolds number must account for the shock moving upstream. 

Reduction of the stagnation pressure also moves the shock up- 
stream as shown in the sequence of photographs in Figure 15. A strong 
shock exists at the nozzle exit of Figure 15a with a stagnation pressure 
of 65.35 inches of mercury. Reduction of the stagnation pressure to 
59.96 inches of mercury allowed the shock to move upstream one and one 
half inches towards the nozzle throat as shown in Figure 15c. Reduction 
of the shock intensity may be noted also as the shock moves upstream by 
comparing Figures 15a and 15c. 

The minimum starting area of the variable diffuser was found to 
range between 0.695 and 0.700 square inches as compared to the isentropic 
minimum area of 0.665 square inches. The area actually required corres- 
ponds to a diffuser heights of 1.39 to 1.40 inches. Figure 16a is a 
photograph of the diffuser nozzle set at a height of 1.39 inches. The 
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shock could not be blown through the diffuser with a stagnation pressure 

of 82.71 inches of mercury. Increasing the diffuser height to 1,40 

inches, (area change of 0.005 square inches) permitted the shock to be 

blown through the diffuser nozzle and the stagnation pressure reduced 

* 

as shown in Figure 16b. The stagnation pressure could have been reduced 
to move the shock farther upstream to the throat position. Figure 17a 
and 17b depict tunnel operation with the diffuser height set slightly 
greater than the required minimum and show the shock standing in the 
diffuser throat and the stagnation pressure reduced to a minimum while 
maintaining the shock location. Several minutes after Figure 17a was 
taken the stagnation pressure was reduced from 60.00 inches of mercury to 
59.50 inches of mercury and the shock moved through the diffuser nozzle 
and relocated itself in the diverging portion of the tunnel nozzle. 

The sharpness of the Schlieren photographs was poor due to the use 
of non-collimated light and the fact the light beam passed through the 
test region twice, Appendix E. Figure 16c gives an indication of the 
seriousness of this problem. The knife edge was not inserted into the 
light beam. The black parallel vertical lines are images produced by 
placing fine copper wires on the outside of the Schlieren window nearest 
to the spherical mirror; the diagonal parallel lines are images produced 
by placing the same size wire on the outside of the Schlieren window 
nearest to the light source. The greater the distance between the 
images of any wire, the farther the light beam had to travel before 
reaching the spherical mirror. The non-parallel light beam intercepted 
the diagonal wire first and hence the image of this wire was deflected 
more due to the longer distance to the spherical mirror than it was for 
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any of the vertical wires which were closer to the spherical mirror. As 
evidenced by the two images formed for each wire, the Schlieren picture 
of the of the air flow would tend to be fuzzy due to superimposing the 
two images formed. The horizontal black line is a wire placed at the 
tunnel centerline. 

Also contributing to poor quality of the Schlieren pictures was 

the incomplete removal of oil vapor from the air. During tunnel opera- 

s 

tion, oil would be deposited in a very thin film on the glass wall and 
cause a refraction of the light. In addition, the film deposited was 
not of a uniform nature. This condition may be seen in Figure 17c. 

The dark convex lines at stations 4 and 5 are the result of an oil film 
deposit. These same lines could still be seen after securing the air 
supply and leaving the Schlieren unit in operation. Due to this oil 
film formation, it was necessary to remove the Schlieren window at 
frequent intervals and wipe off the small amount of oil that had been 
deposited. Water vapor in the atmosphere condensing on the outside 
surface of the Schlieren window was also a factor in obtaining clear 
Schlieren photographs. 

As the vertical wires were spaced at one inch intervals, the dis- 
tance between wires in the photograph may be measured and comparison of 
this measured value to the known one inch dimension permits determining 
a scale factor for the photograph. This method was utilized to deter- 
mine the diffuser heights for heights other than 1.718 inches. A 
lucite gage block, 1.718 inches high, was utilized when runs of this 
diffuser height were made. 
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The apparent difference in intensity of the shock and 
the visibility of the boundary layer caused by the knife 
edge orientation demonstrates why a horizontal knife is 
generally used when studying boundary layers and a vertical 
knife is used when studying shock waves. 

B. Quantitative Results . 

The quantitative evaluation was derived primarily from 
pressure surveys. The Mach numbers actually obtained in 
the tunnel were generally less than the predicted adiabatic 
Mach numbers due to friction effects not accounted for, ob- 
lique shocks and other departures from a one dimensional 
flow analysis. Figure 20 is a plot of the adiabatic Mach 
numbers obtained from the pressure and area ratios for a 
constant diffuser height of 1.718 inches vs. distance down- 
stream of the tunnel nozzle throat. The solid curve is the 
isentropic Mach number based on the change of area ratios 
only as given in Table I, while the dashed curve is the 
predicted adiabatic curve as determined by the use of iso- 
clines in Appendix F. The circles are the adiabatic Mach 
numbers obtained from the data taken on the tunnel and are 
tabulated in Table V, based on the procedure outlined in 
Appendix F. 

Up to a Mach number of 1.936 at pressure tap number 
4, the actual Mach numbers followed the adiabatic curve 
closely, exceeding both the isentropic and adiabatic 
curves between pressure taps numbers 4 and 7. At this 
point, the Mach numbers obtained in the tunnel began to 
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fall below the isentropic curve but still above the adia- 
batic curve to pressure tap number 9. The maximum Mach 
number of 2.275 was obtained at pressure tap number 7. The 
decrease in Mach number at taps number 10 and 24 is un- 
explained. This decrease in Mach number at these pressure 
taps appeared consistently in all runs made during the 
evaluation of the tunnel. 

An attempt was made to maintain a constant Mach number 
in the test section by increasing the height of the varia- 
ble area diffuser. The results of this run are shown 
graphically in Figure 21 and in tabular form in Table VII. 
Because of design limitations, the diffuser height could 
not exceed 1.883 inches. The calculated maximum height of 
the diffuser to maintain a constant Mach number in the test 
section including the effects of friction was 2.145 inches, 
Appendix F, based on a friction factor of 4f = 0.01. Table 
IV presents the diffuser coordinates vs. diffuser length in 
tabular form while Figure 19b presents the same data in 
graphical form. Examination of Figure 21 indicates a 
definite drop in Mach number between the nozzle exit and 
tap 24, in much the same manner as the drop for constant 
diffuser area. Figure 20, between the same pressure tap 
locations. As the diffuser did not increase in height 
appreciably in the first 1.5 inches of diffuser length, the 
effect of friction was to continually reduce the Mach num- 
ber. Recovery of the Mach number was accomplished shortly 
after the diffuser height began to increase. The recovered 
Mach number is somewhat higher than that for constant 
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diffuser area, and appears to be maintaining a more con- 
stant value, thus indicating the effect of increasing the 
area of the diffuser. Had it been possible to increase the 
diffuser height to the calculated height of 2.145 inches, 
the resulting Mach number would have been increased much 
more. Until the behavior in Mach number without area 
change can be explained, the effect of changing the area of 
the diffuser can not be adequately predicted. 

Figure 24 is a plot of adiabatic Mach numbers vs. 
distance from the tunnel nozzle throat for two operating 
conditions: (1) the diffuser height set at constant 1.718 

inches and a back pressure of 38.98 inches of mercury, and 
(2) the diffuser set at the minimum starting height of 1.40 
inches with a back pressure of 30.54 inches of mercury. 

The minimum diffuser height was based on the actual minimum 
height required to operate the tunnel supersonically in the 
test section confirmed by the calculated height in Appendix 
F. The calculated minimum height for isentropic flow was 
found to be 1.33 inches; while that for adiabatic flow was 
1.40 inches. Appendix F. Table III gives the diffuser co- 
ordinates in tabular form and Figure 19b presents the same 
information graphically. It should be noted there was 
very little change of area of the diffuser near the nozzle 
outlet when adjustment of the diffuser height was made 
downstream. 

The data followed the predicted curves in much the 
same general pattern as when the tunnel was operating under 
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a constant diffuser height shock free condition. The 
characteristic drop in Mach numbers may be noted at pres- 
sure taps 10 and 24. A predicted adiabatic Mach number at 
the diffuser throat was calculated with the diffuser height 
set at 1.40 inches, Appendix F. This predicted value of 
M t= 1.64 does not correlate too well with the value of 
M = 1.34 obtained in the tunnel as shown in Figure 24. 

While operating supersonically, the tunnel had a mass 
flow rate ranging from 0.25 Ibm/sec. to 0.36 Ibm/sec., as 
determined by the use of an orifice meter. The correspond- 
ing stagnation pressure varied from 62 to 88 inches of 
mercury. The mass rate of flow for the orifice meter may be 
determined by the procedure outlined in Appendix C, or the 
ideal mass rate of flow for the tunnel nozzle throat, as 
determined by Fliegner's formula, may be corrected by the 
use of the empirical Coefficient of Discharge presented 
graphically in Figure 25. 

The average friction factor 4f, in the diffuser sec- 
tion was found to be 0.01865 when the tunnel was operating 
shock free supersonically. Appendix F. The diffuser height 
was maintained at a constant 1.718 inches and the data 
from Table V utilized. All predictions made involving the 
use of friction were based on the assumption the mean 
friction coefficient 4f was 0.01 and was constant at this 
value. 

Figures 22 and 23, and Table IX present the results 
of a pressure traverse of the tunnel operating with a 
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normal shock in the tunnel nozzle. Figure 22 is a plot of 
the Mach number vs. distance from the nozzle throat. The 
solid line represents the isentropic case assuming a single 
normal shock occurring at M = 2.05 and the dashed line the 
adiabatic case assuming a normal shock at M = 2.00. Imme- 
diately downstream of the single shock the isentropic Mach 
number will be decreased at 0.569 and continue to decrease 
in the nozzle, due to the increasing area, until it reaches 
a value of 0.388 at the end of the nozzle; the isentropic 
Mach number then remains constant in the diffuser section. 
For the adiabatic case, the Mach number immediately down- 
stream of the shock will be 0.577, based on the shock occur- 
ring at M = 2.00. The Mach number continues to decrease for 
the subsonic flow but at a slower rate. Between pressure 
taps 9 and 15, the area change is slight and the Mach number 
remains constant at M - 0.429. In the constant area dif- 
fuser, the adiabatic Mach number is increased due to the 
frictional effects. As shown in the Schlieren photographs 
by the repeated shock pattern. Figures 13 to 18, the flow 
conditions through a normal shock are definitely not a 
simple one dimensional process. A shock occurs and the flow 
immediately goes subsonic. Due to inner action with the 
boundary layer, the flow accelerates to a lower supersonic 
value and the process is repeated until the flow remains 
subsonic. The overall effect of Mach number reduction 
caused by the multiple shocks is shown in Figure 22. Due 
to the complex inner action of the normal shocks and the 
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oblique shocks with the boundary layer, and the change of 
stagnation pressure for each shock occurring, definite 
determination of the Mach number through a shocking system 
is difficult. The method for calculating the Mach numbers 
of Figure 22 is presented in Appendix F. 

At a Mach number exceeding 1.00 the effect of friction 
is to reduce the Mach number to a limiting value of 1.00. 

At Mach numbers lower than unity, the effect of friction is 
reversed; namely, the subsonic flow is accelerated to a 
limiting Mach number of 1„00. This increase of Mach number 
due to the effects of friction on subsonic flow may be noted 
by the gradual increase of Mach number of the subsonic flow 
in Figure 22 or in Table IX. 

Figure 23 is a plot of the tunnel pressure ratios 
based on the stilling tank stagnation pressure vs. distance 
from the nozzle throat for the cases of shock free and 
shock in the nozzle operations. The pressure ratio pattern 
for shock free operation follows the isentropic curve 
closely with deviations commencing about four inches from 
the nozzle throat. With shock in the nozzle operation, 
failure of the pressure ratio downstream of the shock to 
recover to the isentropic value is a loss of stagnation 
pressure caused by the repeated shock pattern. 
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Conclusions and Recommendations 

A. Conclusions . 

1. The qualitative and quantitative evaluation of the tun- 
nel system as installed, show it to be a suitable laboratory 
device for the study of compressible fluid flow. 

2. The minimum starting height of the variable area dif- 
fuser was found to be 1.40 inches. 

3. The average friction factor, 4f, was found to be 0.02044, 
but is believed high in value. 

4. Schlieren photographs of the flow field may be taken in 
black and white or in color. 

5. A vertical knife edge orientation in the Schlieren sys- 
tem is best for studying shock phenomenon. 

6. A horizontal knife edge orientation in the Schlieren 
system is best for studying the boundary layer. 

7. The sharpness of the Schlieren photographs is poor due 
to the use of non-collimated light and the light beam pass- 
ing through the tunnel test section twice. 

8. The phenomenon of a standing shock is definitely not 
one dimensional. 

B. Recommendations . 

It is recommended the following modifications for im- 
proving the performance and usefulness of the tunnel system 
be undertaken : 

1. An additional relay be installed on the Bailey stagna- 
tion pressure controller to maintain a more uniform stag- 
nation pressure. 
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2. Additional heat capacitors be installed in the receiv- 
ers to maintain a more constant stagnation temperature. 

3. A total pressure probe be installed in the fixed geom- 
etry diffuser to obtain the back pressure stagnation pressure. 

4. For more quantitative visual results, the Schlieren 
system be modified to a single pass parallel light beam 
system. 

5. A more rapid means of obtaining the differential pres- 
sure measurements be installed. 

6. The camera system be modified to a reflex type system. 

7. An additional oil fi],ter be installed. 

8. Additional pressure taps be installed downstream of the 
variable area diffuser in the existing pressure survey tun- 
nel sidewall. 

9. The geometry of the variable area diffuser be modified 
to permit constant Mach number runs. 
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Figure 5. Stilling Tank. 




Figure 10. Manometer Manifolds 
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Figure 6. Test Section, Disassembled 
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Figure 7. Stilling Tank, Entrance Nozzle, Test Section 
and Diffuser Assembly. 
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Figure 9. Overall View of Tunnel Section from Front 
Side. 
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Figure 8. Overall View of Tunnel Section from Rear 
Side. 



I 







k) 

>ci- 



k) Uj 
U. 




30 



FI6. II . 5c/iUEREN Camera 5 r3TEM 



Figure 12. Schlleren System Light Source, Knife Edge, 
and Camera -Viewing Components. 





A. VtrCLcal Antfe. P<, - 71.46 In Hg. Pfe • 30.92 to. Hg. - 61.1* F. 

Nozzlt RunnlDg Shock Frca wLCh Atr Flow RlgbC Co LafC. Rota Strong Mach 
Uavalac OrLgLoaclng ac JuncCtnn of Nnzzla and TaaC SacCtnn. Saa Flgura 18. 




B. Vartlcal KnLfa. Pq • 72.31 In Hg. a., Pj, • 30.92 In. Hg. a., • 59.23 P. 

Tonnal Running Shock Frae Wlch Atr Plow Right to Left. 





C. Horizontal Kalfa. P^ - 72.71 In. Hg. a., Pw - 30.92 In. Hg. a., Tq - 65.5* F. 
Tonnal Running Shock Fraa With Air Flow Right to LafC. Kota Thlcknaaa of Cha 
Boundary Layara. 



Tonnal Panning Shock Praa 
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Flgura 13. 





A. Varclcal Knife. 
Dlffuaer Haight 



■ 67,35 In. Hg. a., - 30.45 In. Hg. a., T j 66.3* P. 

1.718 Inchea. Nnrmal Shock Station 12. 




8. Hnrltontal Balte. 
Dlffuaar Haight 



P. • *7.41 In. Hg. a., m 30.86 In. Hg. 
1.718 lachaa. Hotval Shock at Station 12. 



- 65.9' r. 




C. Vartlcal Knlfa. - 81.52 In, Hg. a., - 38.32 In. Hg. a., T . 57.0* P. 

Dlffuaer Haight 1.716 Inchaa. Dua to Higher Danalty and Rsynolda NoBber 
Shock haa Hovad Dpatraam aa Compared to Figure 14A. 



Figure 14. Effect of Incraaaad hack Praaaura on Shock Location 
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A. Vertical Knife. - 65.35 In. Hg. e., P^ - 30.65 In. Hg, e., T_ • 59.6* P. 
Boundary Layer Vlalbla on Nnzzla Walla and Incraaalng In Thlcknaaa Downatraaoi 
of Shnck. 




u *’o “ "*• •** »«• To - b3.2* F. 

Note Pull off of Boundary Layer at Shock and AaynMtrlcal Shading of Cnrraa- 
pnndlng Parta on Oppnalte Sldaa of Tunnel Centerline. Shock la MoTlng Tnvarda 
Noetle Throat aa P„ la Decraaaad fro* Figure 15 A. 




C. Vertical Knife. - 59.96 lo. 
Weaker Shock Moving Opatrea* aa 



Hg. a., P,, - 30.73 In. Hg. a., T„ - 61. g* 
pQ la Raducad an Additional Anwunt. 



F. 



Figure 15. Effect nf Reducing Stagnatloo Praaaura for Shnck lo Nnacla 
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A. Varttcal Knife. 
Shock Could not 



?o • 82.71 In. Hg.e., - 57.6* 
be Blown Through Dtffueer Notsle. 



r, Diffueer Haight 




B. Vertical Knife. Fq - ^0-26 In. Hg. a., - 58.6* F. Olffoeer Haight - 1.40 to. 

Araa haa been Incraaaad 0.005 Square lochaa Conparad to Figure 18 A, aod Shock 
haa bean Blown Through Dlffuaar Nozzle. Praatora Coold be Reduced o»re to H>ve 
Shock towarda the Throat. 




C. Tunnel oot Running. Knife Edge not loaerted. Effect of Noo-Perellel Two Peaa 
Light la Evldaot. Double Vertlcel Lloea era iMgea of Ulraa Placed oo 
Spherical Mirror Side o&.Schlleren Wlodow. Dlegonel Ltoea ere leiagea of Wlrea 
Pieced on Light Source Side of Schlleren Window. Light Peaalng Theae Wlrea 
had to Travel Farther before .Striking Spherlcel Mirror end Hence More Deflection 
of the Tnaee. 



Figure 18. Effect of Aree Change on Shock Location In Dlffuaar Noztle. 
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A. V«rtlc«l Kniffl. Pg - 60.00 la. Hf. > S5.0* F. Dlffasar 

RfllCVit • 1.42 iacbct. Shock Failed Sack Throogh Olffiuor Rossle 
wheo Pa Deere* ted to 59. SO la. Bg. e. 




>. Borleoatel Ihlfe. P^ • 60.61 lo. Hg. *., • 63.1* F. Dlffaeer 

Height • 1.47 laehee. Rote Streag Pall Off of Boaadery Layer et 
Statloa 12.25. 




C. Vertical Kalfe. P. • 72.31 1a. Rg. *., To • 58.7* F. Dark Coacaae 
Llaee at Statlooa 4 ead 5 era the Haealt of Oil Flla Depoalt oo Ualle. 
Dark Berlsoatel Liaee Starting et Station 3 ere elao oil Streak* on 
Sehllarea Halla. 



Figaro 17. Taaaal taaaiag With Radaced Dlffaaer Area. 
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A. Vertical Knife. - 71.46 lo. H». a. - 30.92 lo. Hs* •. T • 61.1* F, 
Air Flow la Xlght to Uft. Section la ao Eolargenant of Flgora 13 A. 

Nostla la Xnnlng Shock Free. Mach Nwbar datamload fro* Wave let a by 
ReUtlooablp SfO ^ .At Corraapoodlng Polot Adiabatic Mach 

Nwbar la 2.125. 




B. Vertical Knife. P^ - 72.31 lo. Hg. Pj, • 30.92 lo. Hg. a., - 59.23* F. 

Air Flow la Right to Left. Section la an Eolargaaiaot of Figure 13B. At 
Corraapoodlng Point Adlebetlc Mech Ntaiber la 1.828. 



Figure 18. Datannlaatloo of Mach Rtaaber froo> Mach Wavalata. 
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A. Vartlc«l Knlf*. - 71.46 to. Hg. «. P^ - 30.92 to. Hg. a. T - 61.1* P. 
Air Plow la Ktghc to Uft. Sactton la ao Bolargaatant of Flgara 13 A. 

Noxala la Rannlog Shock Fraa. Hach RtMibar daCamitnad from Wavalata by 
Ralatlooahlp S//? ^ .At Corratpondlng Point Adiabatic Mach 

tf»bar la 2.125. 




B. Vartlcal lolfa. P„ - 72.31 In. Hg. a., P^ - 30.92 In. Hg. a., - 59.23* P. 

Air Flow la Right to Uft. Sactloo la an Enlargaoant of Flgora 13B. At 
Corraapondlng Point Adiabatic Mach Numbar la 1.628. 



Flgura 18. Datarmlnatlon of Hach Htaabar fron Hach Wavalata. 



9 




3 9 



I 

I 




i 

< 


2.340 

2.337 

2.335 

2.332 

2.332 

2.308 

2.262 

2.215 

2.165 

2.115 

2.090 

1.980 

1.943 

1.920 

1.908 


CM 

<5 


P.sssssSSSISSgl 

OOOOOOOOOOOOOOO 


Distance 

from 

1 exit 




ilsIsllsllllH.I 


H 




/it 


0.0716 

0.0718 

0.0771 

0.0851 

0.0935 

0.099*^ 

0.1021 

0.0956 

0.0899 

0.0850 

0.0820 

0.0762 


M 

C topic 

! 


SI.5?sSiiSssS 

cj CM CM r>i CM <N M r» «M «M ft 


A/A* 


2^340 

2.332 

2.240 

2.U5 

2.005 

1.93? 

1.908 

1.979 

2,050 

2.115 

2.190 

2.255 




sHcSIsSScis 


0OOC3O0OOOO0O 


Distance 

from 

Nozzle 

exit 

in. 


sgssssssssss 



A/A* 


sSS5S3S|5K2§s3a 

^ ^ ~ ~ ^ fj fj rJ rJ fsj* rJ rJ rJ 




< 


ssSssSSisHSIIs 


OOOOOOOOOOOOOOO 


Distance 

from 

exit. In. 


doood-*-'<^fM»nrn>j.jinm 


H 


72SSSSSaSSSSSSS 


tropic 


0.0716 

0.0716 

0.0689 

0.0658 

0.0629 

0.0613 

0.0617 

0.0618 

0.0645 

0.0660 

0.0676 

0.0685 


1 

M 

tropic 


CMfJCMCMCMCMCMC'ICMfMfMCM 


1 

K/K* 


2.340 

2.339 

2.393 

2.460 

2.525 

2.565 

2.555 

2.525 

2.490 

2,455 

2.420 

2.400 


CM 

<i 


sisSssHsiSi 


dddddddooood 


Distance 

nozzle 
exit . in. 


8§§g88g§8§88 





A/ A* 






ooooocooooc:o 


Distance 
throat. In. 

-1.601 

-1.101 

-0.601 

-0.101 

0.399 

0.899 

1.399 

1.899 

2.399 
2.649 

2.899 
3.149 


H 


CM 


s'ij 


dddddddddddddddddddd 


■i! 


, §SS2S»SSI§3l5l5!l°ss 




A/a* 


iissssisSSSSSISSSSSS 






SSSIgSl^ssilsssssss 

o d d d d o d o o d o d o o d d d d d d 


Distance 
throat. in, 


.000 

.171 

.343 

.604 

.686 

.857 

,029 

.200 

.372 

,543 

.715 

.886 

.058 

.229 

.401 

,562 

,743 

,915 

.086 

.201 


O OOOOO-^^-^^-^-^fMCMCMCMCMCMrSr^ 





1.000 

1.097 

1.3m0 

1.621 

1.883 

2.060 

2.190 

2.271 

2.315 

2.330 

2.335 

2.240 


CM 


0.367 

0.402 

0.492 

0.595 

0.691 

0.756 

0.804 

0.834 

0.850 

0.855 

0.857 

0.959 


Distance 

from 

Throat . 

O.COO 
0.500 
1.000 

1.500 
2.000 

2.500 

3.000 

3.500 

4.000 
4.250 

4.500 
4.750 




oo 

^CMt-l.JtrtsOr.OOO'O^^ 

r,j 




iliisiissilliiiilsiili 


dddddodddddddododoodoo 




,000 

,157 

,315 

.478 

,636 

,778 

,894 

,987 

.069 

.133 

,185 

.225 

.264 

.291 

.315 

.331 

.345 

,354 

.364 

.367 

.369 

.271 






1 ^v/v 
1 


1.000 
1.019 
1.073 
1 162 
1.280 
1.415 
1.548 
1.670 
1.788 
1.890 
1.978 
2.050 
2.120 
2.175 
2.220 
2.255 
2.285 
2.305 
2.320 
2.330 
2.337 
2.340 


CM 


sSS2?5s5sl55Sssssis££s 


ddoododdddddddodoooooo 


Distance 

from 

Throat. in 


i~5SI~5sgss5es5si2s5si 


OOodd'^^^^CMCMCMrMCMC-^c'^f-ic^^'t'J'J 



4 0 









I 




<■1 



r 


2.064 

2.179 

2.064 

2.106 

2.042 

2.001 

1.920 

1.862- 

1.979 

2.023 

2.041 

1.945 

1.909 

1.958 

1.975 

1.753 




iliilssiisiisSsi 


oooooooooooooood 




sIsssssrIISSsISS 


< 




o 

Cl. 


0.0880 

0.0813 

0.0880 

0.0855 

0.0895 

0.0922 

0.0976 

0.1042 

0.0909 

0.0868 

0.0855 

0.0895 

0.0908 

0.0869 

0.0856 

0.1018 


cu 5*i 

<3 




S 2 S S ts S S S 5 5 2 S S 


^ i 




X 


.352 

.613 

.916 

.169 

.265 

.227 

.197 

.190 

.011 

,155 

,085 

.126 

.128 

.103 

.126 

.223 




2 


aJ 


SSS2R222SSg22ls2 


doddddoooddododd 


♦ 

< 


sI5sIS§SSS|ss558 


< 




O 

Cl. 


0.3345 

0.2175 

0.1413 

0.1016 

0.0871 

0.0841 

0.0828 

0.0S16 

0.0921 

0.0828 

0.0868 

0.0843 

0.0841 

0.0857 

0.0843 

0.0843 


CU si 




< 


5222222222222222 


c i 




m 


2222222322222222 


J ii 




o ee 
H o 


SSkSSSSkSSssSSSS 


1 i 


90 

131 

171 

217 

261 

287 

319 

367 

397 

434 

465 

494 

530 

559 

589 

677 



z 


sSSssissSsisSsss 




{MCM<NrLjCMrU.-l..u^..u..u*-..u»-^»u 


Q.\0f 


gSSSSSSSSRSSsSIS 


ddooddddooodddoo 


* 

< 


.340 

.340 

.338 

.335 

.333 

.333 

.308 

.262 

.215 

.165 

.115 

.090 

.980 

.943 

.920 

.908 


< 








Cu 


oooooooooooooooo 


S ii 


c'l.Or^ro.-oOxOr^u^oocMO'Jc^^'a 




S i 




X 


1.359 

1.606 

1.835 

2.166 

2.235 

2.198 

2.204 

2.176 

2.056 

2.135 

2.072 

2.099 

2.156 

2.099 

2.099 

2.220 


^iSiT 


SSS5?;Sg§?lS22?i2Sg 

SSS2222SS2SS2222 




ooooddooodoooooo 


♦ 

< 


sl5il§SSsSnis§S 


< 


.-^^^rUfSIC-Urg<NI(NlfMCNICSJ(VIfN(NI 


o 

5: 

Cu 


0.3320 
0.2180 
0.1502 
0.1018 
0.0889 
0.0857 
0.08 22, 
0.0824 
0.0889 
0.0840 
0.0875 
0.0859 
0.0826 
0.0859 
0.0859 
0.0845 


S 




5 5 IS s IS 


B ^ 




2 




Cu° 




^ tt, 


SskssSSSSSSkSSsS 


1 i 


-S25SS3S5SS|g5SS 



z 






0.2020 

0.1800 

0.1992 

0.1962 

0.2025 

0.2125 

0.2225 

0.2415 

0.2230 

0.2160 

0.2260 

0.2420 

0.2350 

0.2450 

0.2545 

0.2960 


♦ 

< 


nntnnnn%n 


< 


fMPJ<Mf>ICMC>JtM<M<NPvlCM<NrMCMCMfM 


o 

£: 


iiiliiiiiliiiiis 


Cu 


ddddoddoddoooooo 


^ a 




SSS552S3SSS33S33 






z 


1.346 

1.608 

1.936 

2.191 

2.272 

2.275 

2.272 

2.213 

2.001 

2.199 

2.124 
2.172 
2.195 
2.103 

2.125 
2.294 




sisisisissmiss 




dddddddddddodddo 


A/A* 


1.097 

1.340 
1.621 
1.883 
2.060 
2.190 
2.270 
2.315 
2.330 
2.335 

2.340 
2.340 
2.340 
2.340 
2.340 
2.190 


o 


IIEliiiiiiliiili 


Cu 


dddooddddddoddoo 


‘ ii 




5 s s 5 5 s ;s s :s 


Bi 




In. 
He. a. 


S S S S S 5 ? 5 s s s s s s s 




skSkSSSSSkSSSsss 


J s 

H 


0 

39 

75 

110 

150 

189 

236 

272 

299 

340 

379 

415 

441 

477 

516 

550 



z 


sSISIsssissSisRs 




iii!i;sSi6ML58| 




doo'dddoooooooooo 






< 


r»r>»«N«N<M«MIV»CM«MC'><N<MCMCN<MfS 




iiiiliilliiiiils 


Cu 


cdoooodbooooocoo 




0'^oc^O'>£>f''ao.-‘c»^-<^ r~0'<oriO' 


^ ii 


S S g g g g s s: s s s s ;e 


^ £ 






,350 

,619 

,933 

,183 

,252 

,254 

.243 

.193 

.998 

.162 

.095 

.113 

.152 

.111 

.150 

.269 




.^.^^rsirurMr'irM.-'Cicucurutsirjru 


:=*-=r 


slSIlliilsilsIll 




ddooodddodoooooo 


< 


097 

340 

621 

883 

060 

190 

270 

315 

330 

335 

340 

340 

340 

340 

340 

190 




.I^rUr^CM*C>JrJrJr'if>Io4cJcMCMfv4CM 


o 


IlSIiiliirilliil 


Cu 


OOOOOOOOOOOOOOOO 


^ is 


— 00>00»Siri>0<*</^M> 


ssggssssgsssssss 


Tap 

No. 


<Nf^>»kr(>or.ooa>o — 




S2SSS222§gS^ggK2 


i i 


SSSSSSSSSSRSRRSS 


o'" 


526 

525 

524 

523 

5^1 

521 

520 

519 

519 

518 

518 

517 

517 

516 

516 

516 


1 i 


"gSIgSssJlRRsliS 



4 / 



For diffuser areas see Table III 








jj|s? 5i!5;IISSsSSssssSsss 


1 1 loooooooooooodddooodooo 






ooooooooodddddddddddddodd 


X 


sHsSssigIssssIsHHHHs 




r^r-iHc<iHHOOodddddddddddddddd 


-tW 


HH^ssSHigsHlsssssssss 


TJ 


odddddddododddd^H^HHcHpHpH^FH 




sssilsSSsSsIsIsSpips^ 

c-trHrHiHOIfNCttNfMOIOldotCNOKNdddciddoidd 


i 

n 1 


2 



ii 



Sc; 

'I 



s' 


0.3720 

0.376S 

0.3860 

0.4020 

0.4130 

0.4210 

0.4300 

0.4375 

0.4390 

0.4400 

0.4420 

0.4420 

0.4400 

0.4390 

0.4375 

0.4375 


‘is 






Vs 




s' 


0.3320 

0.2130 

0.1358 

0.2110 

0.2115 

0.2305 

0.2560 

0.2920 

0.3280 

0.3420 

0.3780 

0.3735 

093670 

0.3680 

0.3680 

0.2315 


‘i 4 






'is 






j 4 


63.28 

63.18 

63.18 

63.18 

63.08 

63.08 

63.08 

63.08 

63.08 

63.08 

63.08 

62.98 

62.98 

62.88 

62.88 

62.88 






2222222222222222 




21 


SliSsisSSSSSSSSI 



i 





12 




1 

2 



I 



APPENDIX A 



TUNNEL SYSTEM OPERATION PROCEDURE 

I. For charging the storage tanks the following procedure should be 
followed; 

1. Open the 3/4 inch valve connecting the trap with the separa- 
tor. 

2. Close the 3/4 inch drain valve located below the tee after 
the trap. 

3. Open the 1/4 inch balance line valve between the separator 
and the trap. 

4. Open the cooling water valve to the after-cooler heat ex- 
changer. 

5. Close both two inch inlet valves located at the entrance to 
the low pressure storage tanks. 

6. Open the valve on the low pressure after-cooler. 

NOTE: Items 5 and 6 are safety precautions to prevent high 

pressure air from entering the low pressure tanks. 

7. Close the drain valve on each high pressure tank. 

8. Close the two inch discharge valve on each high pressure 
tank outlet. 

9. Open the inlet valve to each high pressure tank. 

10. Place the unloader on the compressor in the START position 
(vertical with the index head down.) 

11. Place cover on the compressor oil vent. 

12. Open the air supply valve to the Bailey controller (Pressure 
should read 28 psi). 
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13. Adjust compressor discharge pressure for a pressure slight^ 
ly higher than any pressure indicated in the receivers* This 
is done by turning the adjustment knob on the Bailey Con- 
troller counterclockwise until the red hand indicates the de- 
sired pressure* 

14. Connect the cooling air circulating centrifugal blower in 
the compressor shack* 

15. Throw the roaster power switch located in the main bulkhead 
switch panel for the High Pressure Air Compressor to the ON 
position. 

16. In the compressor shack, throw the overload trip arm to the 
ON position and push the START button* 

17. After the compressor has started running^ place the unloader 
valve in the RUN position (handle horizontal with the index 
head pointing towards the compressor.) 

18* Adjust the compressor discharge pressure regulator valve to 
give the desired discharge pressure* (Maximum 500 psi) 

19* Maintain a check on the storage tank pressure as an automatic 
unloader for the compressor is not provided* 

20* Should an excess amount of air be coming out of the drain line 
from the trap, it is probable that additional water is required 
in the trap. Water may be added when the compressor is secured 
through the union connection of the balance line until it is 
observed coming out of the trap drain* 

21, To secure the compressor place the unloader valve in the START 
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position, push the STOP button and remove the oil vent cover. 
An emergency STOP button is also located with the Bailey con- 
troller. 

22. Close the inlet valves to the storage tanks. 

23. Close the balance line valve. 

24. Close the valve connecting the separator and trap. 

25. Close the air supply valve to the Bailey Controller. 

26. Blow down excess air in the lines from the compressor to the 
storage tank inlet valves by opening the trap drain valve. 

27. Secure the cooling water. 

28. The drain on each high pressure tank and the inner cooler 
should be blown at periodic intervals to remove any accumula- 
tion of oil or water. 
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position, push the STOP button and remove the oil vent cover® 
An emergency STOP button is also located with the Bailey con- 
troller. 

22. Close the inlet valves to the storage tanks. 

23. Close the balance line valve. 

24. Close the valve connecting the separator and trap. 

25. Close the air supply valve to the Bailey Controller. 

26. Blow down excess air in the lines from the compressor co the 
storage tank inlet valves by opening the trap drain valve® 

27® Secure the cooling water, 

28. The drain on each high pressure tank and the inner cooler 

should be blown at periodic intervals to remove any accumula- 
tion of oil or water. 
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II. Blow down procedure 

1. Open the outlet valves on each high pressure tank. 

2. Open the air supply valve to the stagnation pressure Bailey 
Controller (Supply pressure should read 28 psi), 

3. Open the quick opening valve at the outlet of the regulator 
va 1 ve . 

4. Slowly turn the adjusting knob of the Bailey controller 
counter clockwise to allow a gradual buildup of pressure in 
the stilling tank to the desired stagnation pressure. A load- 
ing pressure of at least 5 psi is required before the re- 
gulator valve xtfill -open. 

5. Figure 1-A gives the approximate differential pressure requir- 
ed to obtain the maximum Mach number in the tunnel test sec- 
tion for a given stagnation pressure measured above atmospheric 
pressure in inches of mercury. 

6. The time of each run may be prolonged as well as reducing the 
rate of change of stagnation temperature by allowing the com- 
pressor to continue charging the storage tanks while operat- 
ing the tunnel. 



46 




47 





APPENDIX B 



AIR MOISTURE CONTENT CALCULATIONS. 

Assume that the perfect gas law relationship applies 

PV = mRT 



where P = gas pressure 

V = gas volume f r 

m = lb mass of gas 

R » universal gas constant of gas, 

T = gas temperature, ° Rankine. 

The ratio of mass water vapor in the air to 
is then given by 



(+ - ibf 

Ibm - 

the mass of air 



^ if 

Civ. _ R- T - .Bt 
~ B X P.^ Pa 

/?o T 



where the subscripts ^ refer to air and refers to the 
water vapor. Utilizing the value of R^ = 53.34 and 

R^ = 85.7 

— = 0.<oZZ5 § =0. bZ25 ■ p'' -p 

^ total 

Upon exit from the heat exchanger, the air is saturated at 
500 psia and 60°F. From steam tables at 60°F 

Psat = 0.2563 psia 



and 



n}j 






0.2 Jt 63 _ ._^QQQo^jq !kj2£fS'^ 
S00-O.25h3 ’ Ib.dfyar 
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APPENDIX C 



FLOW RATE CALCULATION PROCEDURE. 

The mass rate of flow for the tunnel, as measured at 
the orifice plate, was determined by the procedure out- 
lined by the American Society of Mechanical Engineers [5] . 
All equations, figures or table numbers correspond to num- 
bers in that publication. Due to the removal of the water 
vapor by the air drying system, it is assumed the air is 
dry for metering purposes. In addition, due to the low 
velocity of the air, the air is considered incompressible. 

For a representative run: 

P^^ = 72.08 in. Hg. 

T = 525°R = 25°F 

o 

Pob = 30.84 in. Hg, 

D = Pipe diameter = 3.120 inches, 
d = 2.496" = orifice diameter 

P orifice diameter = ^ = 0.800 

pipe diameter d 

^P = differential pressure in inches of H 2 O » 4.55 
Equation 2, page 57, gives the mass rate of flow Ibm/hr, 



o 



m 




a 
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Where C 



coefficient of discharge 
F = velocity of approach factor 
Fg = thermal expansion factor 
h^ = differential pressure in inches of H^O 

= specific weight of air at the inlet side of the 
primary element, Ib/ft^ 



From Figure 38, page 67, for a stainless steel orifice plate 
at 65°F 



Fa = 1.00 

From Table 9, page 80, at 65°F 

^ air = 0.075650 Ib/ft^ at 29.92 in. Hg. 

Correcting for increase in pressure 

'f aix = 0.07565 X = 0.0762 Ib/ft^ 

d^ = (2.496 in.)^ - 6.230 in.^ 

Assume a value of K = CF = 0.80 and calculate an approx- 

imate m 



m = 359 X 0.800 x 6.230 

= 1053 Ibm/hr. 



\rT. 
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0.0762 



And from Figure 35, page 62, for 65*^F air = .0000122 

Calculate Reynolds number R from Equation 8. 

Tj _ 0.004244 m _ 0.004244 x 1053 = 117,600 

^ Dg JU 3.120 X 0.0000122 



From Table 5, page 26, or from Figure 1-C of this Appendix, 
for yi3 = 0.800, and R = 117,600 the flow coefficient 

is K = 0.8062 

Corrected calculation of flow rate, therefore is 

m = 359~«-0.8062 x 6.230 * 4.55 x 0,0762 ' 

- 1060 Ibm/hr. = 0.2945 Ibm/sec. 
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The ideal mass flow rate at the nozzle throat is determined 
by Fliegner's formula [2]. 

A*P 



ru; 



Where 



M = 




1.00 




0.532 



^ T ' 

' o 

= ideal flow rate Ibm/sec. 

= tunnel throat area in square inches with 



= stagnation pressure psia 
= stagnation temperature in ®Rankine 



0.532 X 0.367 x 72.08 x 0.491 
->1 525 ^ 



= 0.301 Ibm/sec. 



The nozzle discharge coefficient is defined as the 
ratio of the actual nozzle flow to the flow calculated 
from the isentropic laws for the initial and final pres- 
sures of the actual nozzle. Based on this definition the 
coefficient of discharge is 



m 



D 



m 



actual 

ideal 



0. 2945 
0.301 



0.969 



Figure 1-C presents the Flow Coefficient K plotted 
vs. the Reynolds Number for a flat plate orifice with 
pressure taps at D and 1/2 D for a three inch pipe. This 
figure presents the data of Table 5 (5) graphically for 
values of d/D equal to 0.580 and 0.800, Figure 25 presents 
graphically the nozzle discharge coefficient plotted vs, 
the stilling tank stagnation pressure. Values of were 
determined by the use of Fliegner's formula and by the 
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procedure outlined by the American Society of Mechanical 
Engineers (5) as used in this Appendix. 



52 






I 







I 



I 





f.. 





53 





APPENDIX D 



DESIGN OF DIFFERENTIAL PRESSURE TRANSDUCER. 

Differential pressure transducers were originally to have measured 
the differential pressure between the tunnel pressure and the stagna- 
tion pressure. This method of measurement was discarded due to the 
difficulties of fabricating similar output transducers. 

A successful transducer, Figures 2-D and 3-D, was made based on the 
principle of the deflection of a circular plate held firmly at the edge 
jij . It was desired that each psi of differential pressure produce at 
least 10 microinches per inch of strain. Transforming the equation of 
stress at the center of the clamped plate to the equation of strain (end 
of this appendix) and substituting the appropriate values the plate thick- 

7 

ness was calculated as being 0.0336 inches times the radius of the plate. 
One A-7 strain gage was located at the center of each side of the place 
and the cumulative output of the two gages was read on a Baldwin strain 
indicator. This cumulative output was thus expected to be approximately 
20 microinches per inch of strain per psi of differential pressure. 

The plate thickness actually used was 0.032 inches due to the 
availability of this uniform thickness sheet stock. This plate was 
welded between two hollow cylinderical shells of 1.5 inches inside 
diameter and 2,25 inches outside diameter. IJhile this combination of 
plate thickness and effective plate diameter exceeded the calculated 
plate thickness to plate diameter ratio, the combined output still ex- 
ceeded the desired minimum of 10 microinches. The open ends of the 
cylinderical shells were threaded to receive a plug which contained a 
pressure tap. An 0-ring seal was used between the plug and the shell 
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of the transducer. 



One leg of each strain gage was led through the shell wall and 
the drilled hole secured with epoxy cement. The other leg of the gage 
was soldered to the inside of the shell which served as a common ground 
between the two gages. The two protruding wires from the gages and a 
common ground wire, soldered to the outside of the shell, were then 
connected to the Baldwin strain gage indicator. 

The completed transducer was calibrated on a hydraulic dead weight 
tester. One end of the transducer was vented to the atmosphere and the 
other cavity filled with oil and connected to the dead weight tester. 

A pressure differential ranging from 0 to 100 psi, in increments of 
10 pounds^ was then imposed with strains recorded while loading and un- 
loading the tester. 

A linear plot of pressure differential vs. strain was obtained for 
several transducers but on others, the* welding caused a slight permanent 
buckle in the plate and the strain indications were not reproduceable 
for a given pressure differential. Due to the time involved in making 
and testing each transducer, it was felt that too much unproductive time 
would be spent on making transducers that could not be used; so this 

■ i • 

means of pressure measurement was abandoned in favor of manometer 
pressure measurements. 

One test transducer w^s glided together with Eastman 910 glue to 
eliminate the problems caused by the buckling of the diaphragm due to 
welding. A centering hole and shoulder was made in the cylinderical 
shells, the diaphragm was then glued to both shells under light pressure. 
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The strain gages were attached to each side of the diaphragm prior to 
joining the two end shells. This transducer gave a linear plot of 18 .,52 
microinches per inch of strain/psi of pressure differential when tested 
on the hydraulic dead weight tester but it did not respond to less than 
3/4 of an inch of mercury pressure differential when checked against a 
manometer reading of the same pressure differential with the tunnel in 
operationo 

The stress distribution in a uniformly loaded diaphragm with clamp- 
ed edges is given in Figure 1-D below I'ul . 




Figure 1-D Stress Distributions of Diaphragm 
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The tangential stress at the center of the plate is given by the 



equation 



S = 



Where 



3Ujf£)_Pr 
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S = stress in psi. 

= Poisson ratio, 

|0 = Pressure differential in psi. 

iP = Radius of diaphragm in inches. 

■£ = Diaphragm thickness in inches. 

The corresponding value of strain is obtained from 



( 1 ) 



6 = - - /4 

^ a f E 



( 2 ) 



where ^ is the strain and ^ = modulus of rigidity of diaphragm 
material. Substituting equation (2) into equation (1) 






8 Ef^ 



(3) 



Utilizing the required minimum value of ^ = 10 ^ in^ y' P®i 

differential, ~ E = 30 x 10^ psi for the values of 

steel the ratio of the diaphragm thickness to its radius may be obtained. 



to 






3' r’-j 1-0. 3^ ) 

S ^ 3 k io~i ^ iz. 



i - 0.0 55b 



<4) 



For a diaphragm radius of 0,75 inches a plate thickness of 0.0252 



Inches is required. 
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Figure 2-D. Exterior View of Welded Pressure Transducer. 





Figure 3-D. Interior View of Welded Pressure Transducer. 
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APPENDIX E 



GENEKAL THEORY OF THE SCHLIEREN METHOD OF OPTICAL INVESTIGATION OF 
COMPRESSIBLE FLUID FLOV7 . 

Apart from the conventional methods of experimentally investigating 
flow patterns oy means of pressure and velocity surveys, compressible 
gas flows lend themselves particularly well to optical methods of in- 
vestigation. 

Fundamentally, the optical methods in common use (the interfero- 
meter, the Schlieren, and the shadowgraph) depend on one of the two 
physical phenomena: (1) the speed of light depends on the index of 

refraction of the medium through which it passes, and the index of re- 
fraction of a gas in turn depends upon its density; and, as a conse- 
quence of this first phenomenon, (2) light passing through a density 
gradient in a gas (and therefore through a gradient in index of re- 
fraction) is deflected in the same manner as though it were passing 
through a prism. In a high-speed gas flow the density changes are 
sufficiently large to make these phenomena sizeable enough for optical 
observation. 

The interferometer, based on phenomenon (1), measures directly 
changes of density, and is primarily suited for quantitative determina- 
tion of the density field. 

The Schlieren method, based on phenomenon (2), measures density 
gradients. Although it is theoretically adaptable to quantitative use, 
it is inferior to the interferometer in this respect, and its greatest 
utility is in giving an easily interpretable picture of the flow field 
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together with a rough picture of the density variations in the flow. 

The shadowgraph method, also based on phenomenon (2), measures 
the second derivative of the density (i.e., the first derivative of 

the density gradient). Therefore it makes visible only those parts 

% ■ ' 

of the flow where the density, gradients change very rapidly, and it 

has found its greatest utility in the study of shock waves. 

Of the three methods mentioned, the interferometer yields the 

most Information, and the shadowgraph the least. On the other hand, 

% 

the interferometer is the most costly and the most difficult to operate^ 
whereas the shadowgraph is the least costly and the easiest to operate^ 
Each method therefore has its own useful niche in experimental workj, 
and the choice of method depends on the nature of the investigation. 
Figure 1-E below shows one of the numerous type of Schlieren 
(striae) arrangements. Light from a. uniformly illuminated line source 



rM 




Figure 1-E Schematic Schlieren System 
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of small but finite width is collimated by the first lens and then 
passes through the test, section. It is then brought to a focus by 
the second lens and projected on the screen. At the focal point, 
where there exists an image of the source, there is introduced a knife- 
edge which cuts off part of the light. With no flow in the test sec- 
tion the knife-edge is usually adjusted so as to Intercept about half 
the light, and the screen is uniformly illuminated by the portion of 
the light escaping the knife-edge. When the flow is established in 
the test section assumed ^ for simplicity, to be two imensional, with 
i-.ach light ray passing through a path of constant air density any light 
lay passing through a region in which there is a density gradient nor- 
mal to the light direction will be deflected as though it has passed 
through a prism. Therefore, depending on the orientation of the knife- 
edge with respect to the density gradient, and on the sign of the den- 
sity gradient, more or less of the light passing through each part of 
the test section will escape the knife-edge and Illuminate the screen. 

Thus the S:hlleren system makes density gradients visible in terms 
of intensity of illumination. A photographic plate at the viewing 
screen records density gradients in the test section as different shades 
of gray. 

For the Schlieren figure drawn it is assumed that the flow in the 
test section is parallel and in the x, y-plane, and that the light 
passes through the test section in the z-direction. 

Since the speed of a wave front of light varies inversely with 
the index of refraction n of the medium through which the light travels, 
it follows that a given wave front will rotate as it passes through a 
gradient in n. Accordingly, the normal to the wave front will follow 
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a curved path; since this normal is what we mean by the 
light ray, the latter is refracted as it passes through the 
density gradient. Noting that n is nearly unity for gases, 
it may be shown that 

j 

- q ra d lenf n 
R ^ 

where R is the radius of curvature of the light ray. The 
total angular deflection ^ of the ray in passing through 
the test section of width L is therefore given by 

^ ^ ~ L grad n 

Resolving this into Cartesian components, and connecting 
the index of refraction with the gas density through the 
empirical Gladstone-Dale equation which states that 



<? ' G-D 

where K_ the Gladstone-Dale constant, is constant for a 
given gas, we get the angulerr deflections of the light in 
the X- and y-directions 




L A 



•a 9 

G -‘0 







If the knife-edge is aligned in, the y-direction (i.e. , 
normal to the flow) as in Figure 1-E, only deflection in 
the x-direction will influence the light passing the knife- 
edge. Hence density gradients in the x-direction will be 
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made visible, but gradients in the y-direction will not be 
visible. Similarly, if the knife-edge is parallel to x, 
only gradients in the y-direction will become visible. 

In interpreting a Schlieren photograph of a two dimen- 
sional flow, it is convenient to imagine that the photograph 
represents a view of a relief map of the flow field in which 
the vertical elevation is the gas density. If the map is 
illuminated by glancing side light nearly parallel to x, the 
pattern of shadows and brightness will correspond to a 
Schlieren photograph with the knife-edge in the y-direction. 
Similarly, illumination in the y-direction corresponds to 
placing the Schlieren knife-edge parallel to the x-direction. 

Figures 14b and 14c are Schlieren photographs of the 
installed test section and show the two different effects 
obtainable with orientations of the knife-edge. In Figure 
14b (knife-edge horizontal) light areas on the upper half 
correspond to dark areas on the lower half, and vice versa; 
the boundary layers on the walls are clearly visible, since 
with this knife-edge orientation the density gradients nor- 
mal to the flow are made visible. In Figure 14c (knife edge 
vertical) the upper and lower halves look alike; the bound- 
ary layers are not visible, except aft of the shock wave 
which is inneracting with the boundary layer. 

The quality of focus of the Schlieren pictures is be- 
lieved to be poor due in part to the present arrangement of 
the Schlieren system installed for the reasons as outlined 
below by Holder and North [7] : 

Although it is usually desirable to arrange for the 
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beam of light passing through the working section to be 
parallel, nonparallel beams are occasionally used and an 
arrangement for producing this type is sketched in Figure 
2-Eo Only one mirror or lens is required (apart from the 
focussing lens) , but this advantage is to some extent off- 
set because its diameter must be larger in relation to the 
field to be examined than if a parallel beam is used. 



In addition to employing a nonparallel light beam, 
the mirror system involves a double passage of the beam 
across the working section. The source is placed at the 
center of curvature of the mirror , and part of the image 
of the source is reflected on to the viewing screen by a 
small plane mirror held near the source. Because the 




Figure 2-E. Coincidence Schlieren System using 
a Plane Mirror or knife-edge 
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source and its image are very close to each other, the 
arrangement is often referred to as the "coincidence" system. 
If the image of the source is deflected in a direction per- 
pendicular to the edge of the plane mirror, the illumination 
on the viewing screen changes; the edge of the mirror thus 
plays the same part as the knife edge in a Toepler appa- 
ratus. 

Since the source and its image are not truly coincident, 
the paths of the incident and reflected rays across the 
working section differ even in the absence of deflections 
produced by density gradients. This has an adverse effect 
on the image quality (in addition to any effect associated 
with the use of non-parallel light) since two slight offset 
images are produced. (See Figure 16c) 

These images may be made to coincide at the expense of 
considerably reduced illumination by using the arrangement 
sketched in Figure 3-E. 




Figure 3-E. Coincidence Schlieren System using 
Graded Filter 
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Here a graded filter is placed at 45 degrees to the axis at 
the center of curvature of the mirror „ A condenser lens is 
used to produce an image of the source on the surface of 
the filter, and part of the light from this image is re- 
flected into the mirror » The image of the source formed by 
reflection from the mirror is made to coincide with the 
original image on the surface of the filter, and the light 
transmitted by the filter is arranged to fall on the vaew- 
ing screen in the usual way. The source and xts image nay 
thus be made exactly coincident so that, in the absence of 
density gradients, the incident and reflected beams coin- 
cide. This will no longer be the case when the rays are 
deflected by density gradients, and some loss of im.age 
sharpness, or difficulty in interpreting the observations, 
may then result. 

A further difficulty arises because only one of the 
two images formed by the incident and the reflected light 
can be in focus. It is usually best to focus the image 
associated with the reflected light since the incident rays 
are subject to further refraction during their second pass- 
age across the wind tunnel. 

The principal advantage of the system is that the 
double passage of the light beam, and the fact that the 
source image is formed at the center of curvature of the 
mirror instead of its focus, give increased sensitivity. 

For example, in a two-dimensional flow containing a 
region of density gradient of sufficient extent for the 
incident and reflected ray to remain within it, the 
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angular deflection of the ray will be doubled by the use of 
such a system. Moreover, the displacement of the image of 
the source for a given angular deflection will be twice as 
large when the image is at the center of curvature as when 
it is at the focus. Coincidence systems may thus be val- 
uable where extreme sensitivity is required (for example, 
tunnels working with very low density) , although for or- 
dinary use they are undesirable for the reasons mentioned 
above . 

Part of the increased sensitivity associated with a 
double-passage system may be retained, without the diffi- 
culties associated with the use of non parallel light by 
means of the arrangement sketched in Figure 4-E. Here a 
parallel beam of light is produced by a lens or mirror and 
after crossing the working section it is reflected back by 
means of a plane mirror. 




Figure 4-E. Schlieren System with Parallel Light Beams 
Crossing Tunnel Twice with Increased Sensi- 
tivity. 
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Tinsley Laboratories Inc« furnished the basic components for 
the coincidence Schlieren system as well as the design as 
shown in Figure 11. One unit was the six inch spherical 
mirror with a focal length of 80 inches mounted on a small 
pedestal. The other unit, mounted on a wooden frame with 
leveling jacks, consisted of the light source, knife edge, 
surface mirror and encased in a light tight wooden box the 
camera lens, folding surface mirrors and a six inch ground 
glass view plate. These were mounted on a portable stand 
as shown in Figure 9. 

An Ilex Optical Company Number 3 Acme Synchro Shutter 
was installed between the knife edge and the camera lens by 
means of a machined adapter fitted over the lens mounting 
with the shutter bolted to it. This arrangement allowed 
all of the light passing over the knife edge to enter the 
camera lens system but otherwise maintained the light tight 
integrity of the system. 

The six inch ground glass view plate was removed and 
the back of an Ilex Optical Company Korona camera consist- 
ing of a four inch by five inch ground glass view plate 
and film pack holder fitted in its place. 

The procedure to align the Schlieren system is as 
follows referring to Figure 11: 

1. Locate the spherical mirror as close to the test 
section window as practical. 

2. Locate camera, light source, and knife-edge 
assembly such that the light source lies on an axis 
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perpendicular to the window. The light source should be 
approximately 80 inches from the spherical mirror surface,, 
This may be done by moving the spherical mirror in the slots 
provided on the stand. 

3. With the light source "turned-On" adjust the 
spherical mirror to image the light source at the knife- 
edge. If the image is not in focus with the knife-edge 
assembly at approximately the middle of its slide travel 
move the spherical mirror forward or back to correct the 
image focus. (A small black card is helpful in locating 
images and the light beam directions. ) 

4. Adjust the knife-edge screw such that the "Knife” 
does not intercept the light source image. The light should 
then pass through the camera lens, and an image of the test 
region should appear on the view screen. If the view screen 
image is not centered, adjust the small 1st surface mirror 
and the spherical mirror such that the light passes through 
the center of the camera lens and aperture of the knife- 
edge assembly. If the view screen image is yet not centered, 
adjustment of the three small screw heads at the back of the 
camera will provide the final centering adjustment, 

5. The knife-edge can now be freely adjusted to provide 
the required Schlieren sensitivity. 

NOTE: All mirrors are front surface aluminized and 

should not be touched on the reflecting surface. 
They should be cleaned by dusting with a qual- 
ity camels hair brush. 
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The best photographs of the flow field were obtained 
using the following camera settings: 

1. Shutter speed set at 1/200 of a second. 

2. Lens opening at f-11. This value has nothing to 
do with the actual exposure time of the film. 

This value of the lens opening allows all of the 
light beam to enter the camera and helps reduce 
stray light entering the lens system and fogging 
the film. 

3. Royal Pan film with a tungsten film speed rating 
of 160 was used. 

4. Developing time of eight minutes in total darkness 

in a 1 : 1. dilute mixture of Eastman DK 50 devel- 

o 

oper and water at 68 F. 

Colored Schlieren views of the flow field may be ob- 
tained by replacing the knife-edge assembly with the color 
grid purchased from Tinsley Laboratories. Colored pictures 
were not attempted as a means of printing and developing 
the color film was not available. 

The light source may be either a General Electric 
Company bulb number 1968 or number 1924, or a Westinghouse 
bulb number 1633. All three bulbs produce a line source 
of light. If the General Electric bulbs are used, a dif- 
ferent mount must be made for the bulb. 

Both the light source and the knife edge m.ay be rotated 
from the horizontal position to the vertical position to 
obtain the desired direction of the density gradient of the 
flow field. 
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For the reasons already outlined, the focus of this 
Schlieren system is not as sharp as could be obtained with 
a two mirror one pass light beam system, 

A metal frame was fitted over one of the Schlieren 
windows and small copper wires were soldered into position 
at one inch intervals, starting at the throat position, to 
give horizontal dimensions. 
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APPENDIX F 



SAMPLE CALCULATIONS FOR EXPERIMENTAL DATA AND THEORETICAL 
PREDICTIONS. 

PREDICTION OF ADIABATIC MACH NUMBERS FOR NOZZLE AND CON- 
STANT AREA TEST SECTION OPERATING WITHOUT SHOCK. 

Flow is assumed to be isentropic to the nozzle throat 
and adiabatic thereafter. For ease of calculations, the 
nozzle contour as given in Figure 19a was assumed to con- 
sist of five straight line segments as listed in Table X. 

Utilizing the influence coefficients of table 8.2 (2) 
multiplied by M and designating the product F as outlined 
by Shapiro (2) for the combined effects of area change and 
friction results in the following equation 

dM^ = F -I- F^4f ^ • 

A A ^ f D 

The technique of isoclines as given in Chapter 8 of Shapiro 
(2) was the means of arriving at a solution to this differ- 
ential equation. 

The parameters of the area A and the hydraulic dia- 
meter D may be expressed in terms of the tunnel dimensions. 
Area of the tunnel at any point is the product of the tunnel 
width w and the tunnel height h at the point in question. 

As the tunnel is of constant width, the area at any point is 
proportional to the height at that point and the derivative 
of the area proportional to the derivative of the tunnel 
height. The hydraulic diameter is defined as 4 times the 
area divided by the wetted perimeter and expressed in 
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tunnel dimensions is 



4 w h 
2 (w+ h) 



D = 4 X area 

wetted perimeter 



2 w h 
w+ H” 



As the tunnel width w is constant throughout the length of 
the tunnel at 0.500 inches the hydraulic diameter may be 
expressed 



D 



o(0.500) h 

o:'soo+ h 



2 h 
lt2h 



As assumed, the nozzle segment heights will vary lin- 
early in accordance with the equation for a straight line 
h = mx + h^ where m is the slope of the nozzle segment 
at some distance x downstream of the nozzle throat and the 
nozzle throat having an initial height of h^. 



at X 



h = mx+ h 

' o 

dh= mdx ; dx = — 

m 

L, h = hj^ so the slope m 




*- L 



Combining these relationships and substituting them into 
the original differential equation 



dM^ = 



A ^ 2 



4. -F ^ L C W -f h ; 



w h Ch^ - h^) 



dh 



dM^ = 



dM 



F F^(4f) L 

A , f 

E~ + 



2 w (hi - h ) 



dh 



F F.(4f)L 

-A + _£ 

h 2 w (h^ - h 



^ ) dh 
h 



d+") 



The equation is now in the form so that by suitable numer- 



ical values of h and M being substituted ,^^he slope of the 
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Mach number squared is obtained. Values of are assumed 

in addition to values of M that are thought to bracket the 

2 

actual Mach number. On a plot of M vs, h^ the calculated 

2 

values of the slope dM /dh are plotted at the coordinate 

points of each assumed value of M and h. Sufficiently small 

increments of h must be utilized to provide the necessary 

isoclinic field, A smooth curve , starting tangent to be 

2 

the value of dM /dh at the known value of h^ and its corres- 

o 

ponding value of M is then drawn through the field of iso- 
clines tangent to a reference isocline at each value of h 
the curve passes through. In addition to the physical data 
listed in the previous table a value of the average friction 
factor (4f) equaling 0,01 was assumed with an initial Mach 

number of 1,00 existing in the nozzle throat. 

2 

A value of M is determined from the junction of the 
smooth curve and the height at the segment end point. This 
calculated Mach number is then used as the initial Mach 
number for the second nozzle segment. Continuation of this 
process resulted in the adiabatic Mach number at the noz- 
zle exit. Sample calculations for segment 2 are given in 
tabular form in Table XII and are plotted on figure 20. 

Straightforward use of Table 42(3) for Fanno line 
procedure was used to predict the resultant Mach number at 
the end of the constant area section of the tunnel. Iso- 
clinic procedures could have been used but this method 
was discarded as being time consuming compared to the use 
of the tables. The following information was assurcied or 
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utilized for determining the exit Mach number; a i f 



T 



4f 



I- 



D 

Ml 



0,01 and constant in value throughout the 
section 

5.75 inches, the tunnel section length 
0.775 inches, the hydraulic diameter 
Mach number at the nozzle exit equal 2,17 
Mach number at the section exit 



Entering Table 42(3) with Mt of 2.17, a value of 4fL _ /D 

IHaX 

was read as 0.35298 for use in the following equation: 




0.35298 - (0.01) (5,75) _ 0.27878 

0.775 



Table 42(3) was then entered using this value of 4fL/D and 
the Mach number at the end of the test section determined 
as 1.914. 



PREDICTION OF ADIABATIC MACH NUMBER WITH A NORMAL SHOCK 
STANDING IN THE TUNNEL NOZZLE. 

Flow was assumed to be adiabatic in the nozzle until 
the normal shock occurred at a point midway between pres- 
sure taps 4 and 5. M^ the Mach number upstream of the 
normal shock was determined as 2,0 as read from figure 20, 
which is a plot of the adiabatic Mach numbers as calculated 
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vs. tunnel length based on previous calculations in this 
appendix. Entering Table 48(3) with this value of M , the 
Mach number M^, downstream of the normal shock was obtained 
as 0.57735. thus obtained was used as the initial value 

of M for the isocline technique for determining Mach numbers 
for points downstream of the normal shock in the nozzle and 
in the constant area section of the tunnel. The nozzle area 
and the corresponding nozzle height were read from figure 
19a at the shock location assumed to be between pressure 
taps 4 and 5. Nozzle height at the assumed shock location 
is 1.30 inches. 

As before, the remaining portion downstream of the shock 
in the diverging nozzle was assumed to consist of the same 
straight line segments. The initial section was taken as 
the remaining portion of the original segment 2 and was 
designated as segment 2a. All previous assumptions con- 
cerning the average friction factor and the physical data 
remained constant as outlined in Table XI. 

A nozzle exit Mach number was found to be 0,429. The 
relationship of 4fL/D at the inlet and exit of the constant 
area section was utilized to determine the Mach number at 
the end of the constant area section as was done in pre- 
vious calculations in this appendix. This exit Mach number 
was found to be M = 0,434. 

Theoretical values of P/P , at any point downstream of 

ol 

the shock are obtained by multiplying the value of P/P^ 
corresponding to the Mach number at that point by the value 
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of P /P obtained from Table 48 (3) for the Mach number 

at which the shock occurred. Based on the data in Table 

IX, it was assumed the shock occurred at a Mach number of 

2.000 with P /P having a value of 0.7208. Adiabatic 
oy ox 

Mach numbers and pressure ratios obtained by this procedure 
are presented graphically in Figures 22 and 23. 

Table XVI is a summary of the adiabatic Mach numbers 
predicted for the tunnel operating with and without a nor- 
mal shock standing in the tunnel nozzle with a constant 
area section. 

TABLE XVI - Summary of Adiabatic Mach Numbers 
in Nozzle 

Pressure tap Adiabatic Mach Adiabatic Mach 

Number Number 

No Shock Shock in Nozzle 



2 


1.355 


1.355 


5 


2.070 


0.529 


7 


2.180 


0.432 


9 


2.195 


0.429 


12 


2.170 


0.429 


32 


1.914 


0.434 



(end of constant 
area section) 
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DETERMINATION OF DIFFUSER HEIGHT TO PRODUCE CONSTANT MACH 
NUMBER IN TEST SECTION INCLUDING EFFECT OF FRICTION. 

Assume that the diffuser height will increase linearly. 
Point t is taken as the nozzle throat; point 2, the end of 
the nozzle; and point 3 the maximum diffuser height loca- 
tion. 

Utilizing the influence coefficients for constant 
specific heat and molecular weight from Shapiro's Table 
8.2 [ 2 ] 

' 7 ^'^ "" i-M^ ' A ^ 



Area of tunnel - width x height = w h and as tunnel is of 
constant width 



dA ^ dh Qnd ^ 




h 



The hydraulic diameter D = 



4 (area) 

Wetted perimeter 



_ 4 w h 

2 (W-f h) 



2 w h . Substituting the values 
w + h 



^ A ^ 
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Combining terms 



4- ujdh 
(JJ+h 



= 4 -fJk M^&a 




Integrating between designated limits 



h — 

^2 

In { ^ ^ ^ 



clearing and solving for 

= — (jU -h { LL>+ h ) & ^ 



Substituting values of w = 0.5 inches; Y 12 - 1.718 inch; 

M 2 = 2.17; - X 2 = 5.33 inches. Figure 20. 

Xt(^. i7fia o025)io3 X 

h3 = 0.5 -t (2.218) 0 

h 3 = 0.5 f 2. 645 = 2.145 inches. 

The maximum theoretical height of the diffuser to ob- 
tain a constant adiabatic Mach number of 2.17 could not be 
obtained in the diffuser section due to physical limita- 
tions of the tunnel itself. As may be noted in Figure 19b, 
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the diffuser contours actually existing are far from the 
theoretical solution of the maximum diffuser height. The 
actual contour required is of an exponential nature as 
developed previously. Values required in the diffuser to 
produce a constant adiabatic Mach number in the diffuser 
vs. diffuser length are presented in Table XIII in tabular 
form and plotted in Figure 19b. Numerical values of x were 
inserted in the exponential equation derived for determin- 
ing the diffuser height. 

* 

DETERMINATION OF MINIMUM DIFFUSER AREA TO ALLOW SUPERSONIC 
FLOW, 

As in the case of the diverging nozzle, it was assumed 
the diffuser converged linearly from the maximum nozzle exit 
height of 1.718 inches to some minimum height less than 
1.718 inches. The same isoclinic technique was employed to 
determine this minimum height and allow supersonic flow to 
exist. The most severe condition that may be encountered 
is to have a normal shocTc at the nozzle exit or at some 
other location upstream of the diffuser throat. This con- 
dition will exist upon the initial start of the tunnel blow. 
The minimum diffuser height must be sufficiently large to 
pass the given mass rate of flow in the subsonic condition. 
Assuming a normal shock occurring at the nozzle exit, with 
a Mach number of 2.17 existing in the nozzle exit, the 
initial Mach number for the diffuser will be 0.551 as taken 
from the Gas Tables (3) for normal isentropic shock. At 
the minimum area, the Mach number must increase to 1.00. 
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various values of the minimum diffuser height h are assumed ^ 
the isoclinic field drawn for these conditions and the 
smooth curve faired through the isoclinic field. Various 
values of h, the minimum diffuser height, had to be assumed 
to meet the end point conditions, namely, the curve had to 
have an infinite slope at the junction of h and the Mach 
number squared equaling one. Table XIV presents the data 
for determining the minimum height in tabular form and 
Figure 2-F is a plot of the isdclinic field for a minimum 
diffuser height of 1,40 inches. All isoclines are not 
plotted to reduce the clutter near the terminal point of the 
curves. The slope of dM /dh changes radically near the ter- 
minal point requiring a large number of isoclines to be 
plotted, 

DETERMINATION OF MACH NUMBER AT DIFFUSER THROAT WITH DIF- 
FUSER HEIGHT = 1,40 INCHES CONSIDERING FRICTION 

Conditions at nozzle throat are designated by t, at 
the nozzle exit by 2, and at the diffuser throat by 3, 

Utilizing the influence coefficients of Table 8,2 
2 

multiplied by M and designating the product F as outlined 
by Shapiro [2] for the combined effects of area change and 
friction 




Area of tunnel = width x height = wh 
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As A ^ h for constant width tunnel 





R 4-f dL 




hydraulic diameter D s 


4 X area _ 


4 wh = 


Wetted perimeter 


2 (w-<‘’ h) 


7 wh 






wf h 







w s h inch 

h _ 2 h 

'/e -+h / y- 



Substituting this value’ of 










Assume that the diffuser height will vary linearly in 
accordance with the relationship 

h ~ my. -f Hq 



where m is the slope of the diffuser and x is the distance 
down the tunnel. 

At 



X = = 4,802 inches, h = h 2 = 1,718 inches, 

and at 

X = 10,8 inches, h = h^ = 1,40 inches 




dh mdx 



■= 



_ t.4 7/8 



- - 0. 0530 
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Let 



h = 






and 



Pi - ^^3 



for determination of F 



z 



f - Q ocaa 



/% - /7 from Figure 20 



1 ^ 3 - = R 'v* + 



4f- /> / / 

m 



irn j 






/.55^ 



r , 4-lD.OOzs), 1 

-0.03-30 is.i/e"''/ 



^3 " '''^2 = ^ +G07^^Ff 



Assuming a value of to determine M and the values of 
the influence coefficients from Table B-6 [2] , a process 
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of iteration may be made to find the correct value of 

^ ^ • ' ryi L(o3'i' n _ i 

Trtj A73=r/.63 3 %9 

Fp^ - 4 JKo 55 

^ - /2.038 

- - 0,Z2(F^ *f 

- 4.10- b.e2l(4.1i>35) - OolHa (/2.058) 

3 

A// = ^^70 

A/^ = /.(o4Z 

Trcj 4 M - /• 

Ff - - I^.IO(o 

- ^.70 - 0.22! (4.1019) - 0. 0j92(l2. (OO) 

3 

Ml - 2^9 

-- Lb4 
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In addition to the analytical solution for determining 
the Mach number at the diffuser throat, an isoclinic field 
plot was drawn and an exit Mach number obtained by this 
technique. The Mach number as determined by the use of 
isoclines was 1.71 compared to a Mach number of 1.64 as 
determined by the analytical solution. Table XV and Figure 
3-F present the information necessary to plot the isoclinic 
field and a plot of the field itself. 

DETERMINATION OF EXPERIMENTAL ADIABATIC MACH NUMBERS FOR 
SHOCK FREE CONDITIONS 

Assume that flow to the nozzle throat is isentropic and 
that the flow for the remainder of the tunnel system is 
adiabatic. From the Gas Tables for isentropic flow (3) the 
adiabatic Mach number, M, at any desired point may be deter- 
mined for shock free flow conditions under the product 




providing A* and are taken as the area and stagnation 
pressure when the Mach number is unity. P/Pq = 

1 - — with pressures P^ and AP being measured by the 

use of manometers. Consider pressure tap 16 at the outlet 
of the nozzle as station 1 and tap 31 in the diffuser as 
station 2. 

Data is taken from Table V or Figure 20 for a shock 
free run with the diffuser set at a constant height of 
1.718 inches. 
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2.34 



(va;;) = 

(P/Po)^ = 0.0844 

(A/A^) (P/Pq^i = 2.34 X .0844 = 0.1975 

Prom the Gas Tables for isentropic flow (3) using the above 
value for entering the tables is determined as being 
2.125. All adiabatic Mach numbers upstream of a shock were 
determined in this manner. 



DETERMINATION OF EXPERIMENTAL ADIABATIC MACH NUMBERS FOR 
SHOCK IN NOZZLE OPERATION. 



The adiabatic Mach number for points downstream of the 
shock in the nozzle are determined in the same manner as 
for shock free conditions. Taking data from Table IX for 
pressure tap number 20, the product (A/A*) (P/P^)^ was cal- 

culated as being 0.941. Entering Table 30 (3) with this 
value, the adiabatic Mach number was determined as being 
0.594. The ratio of change in stagnation pressure at pres- 
sure tap 20 or at any other location downstream of the shock 
in the tunnel to the initial stagnation pressure is obtained 
by the relation 



fe). 




n 



The value of (P/Poy^n evaluated from table 30 (3) at 

the Mach number M^ under the column P/Pq- 

DETERMINATION OF EXPERIMENTAL MEAN FRICTION COEFFICIENT 

To eliminate as many variables as possible, data with 
the diffuser set at a constant height of 1.718 inches and 
the tunnel running shock free was used, i.e. , Fannb Plow. 
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With the value of M having been determined from data in Table 
V as being 2.125, at point 1, the constant area inlet, this 
Mach number was used as the argument to enter the Gas Tables 
(3) for Fanno Flow and evaluate the following: 



( P| a 97.7 

^ 0.3407 



Now since 



p ■ {P/P1, 



we obtain 



P 

p*j 



0. 844 ' 



= 0.481 



From the Gas Tables (3) for the Fanno Line flow using 
(p/p*)^ - 0.481 the Mach number at the constant area 

pressure tap 31, point 2, is determined as being 1.78 and 
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the corresponding value of 




0. 23520 



' 4f ^L\ _ 3737 — 0.^350 - 0. 1585 

' D ' 

Q - 0.175 incn&s 



aL = /O.OO" - 4.75" ^ 3.25 inches, 




0J358 « 0.175 

3 . 55 

0, O] S><s5 
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APPENDIX G 



AIR RECEIVER CAPACITANCE DESIGN CALCULATIONS AND PREDICTED 
TEMPERATURE VARIATIONS 

To maintain an isothermal blowdown, Muto 115] states 

if- A 7 cr 

the ratio (I), (7 ~ ^ ^ should have as a minimum 

value of 7.5, and the non dimensional number of transfer 

/? A 



units (2)>NTU = 






should have a minimum value of 7.5, 



Storage tank dimensions 

48 in Outside Diam 

♦ 

11 ft, 4 in^ height 
. 75 in. heads 
. 781 in. shells 
Mass in tanks = Ibm. 



ejm] 

_ r 77 - o.ie/^ / 3 c -f 

■'^WiL me 



- /0j4/0 lb mass of siee-l 



Surface Area of tanks 



-?[7rd . , / g rr = 

u snside a J 



TT(4-(h.44f - 7 

14-4 






fi\REA - 183.5 SCjUar^ foef 



93 



Mass of air initially in tank as beginning of blowdown = Mb 



Initial tank pressure = 475 psia. 



^ pv ^ M4 '%i^ K g ianKi < m 

53, 34 £i -fi X 520'" 

I bm - 

Mb = 577 Ibm. air initially in tank. 

Mass of air in tanks at e(nd of blowdown at 150 psia 

Me . 5 -^1 = l&Z !b„ @ 150 ps/a 

415 ' 

„* . !^le , 182/577 - 0.6! 

^6 

Substituting the above values into equations 1 and 2 



' ^ STBBL _ /Qj 410 K 0 , U 5 „ If n 

^ ' Mi, " 577^ a/75 “ * 

3 ^ ^ HOO Ibm /hr 

h _ 3 X IS3.5 

th 0,173^1100 

2.89 

Since the value of ^ 7.5 and the value of NTU ^7.5 

the value of NTU must be increased by increasing the area 
available for the heat transfer 



using a value of h = 



NTU = 



NTU = 



NTU = 
A = 



7.5 



h A 



~7/r X o. ^ HOP 

3 
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.-1^ »• • i m*4m: 






4T 9^ • » ♦ • 







t^i ^ ^ 




r«l «» «| UJUAi^fp 




I#! 



H 




additional area required 

47b ff^ - fyf 

if 1/8 inch aluminum sheet is used to reduce to corrosion 
problem, the mass of aluminum required utilizing both sides 
of sheet for heat transfer 

„,33 ,,1 . = /b? I? [l.5Z3.ff] 

Mass A1 = 257 Ibm. 

■s 

utilizing eqn 9A as developed by Muto [15] to predict the 
storage tank air temperature after a blow down with the 
presently installed tanks without additional corrective 
heating surfaces and comparing this value to the tempera- 
ture actually noted after a blowdown very close correla- 
tion was noted. 

Eqn 9A by Muto to predict temp after blow 

r / C* \-^ _ 1 0.314 //.S ) 

~ \ !+ Co / ' \ i + /7.8 J 

T* = 0.9785 

with To = 523°R from Figure 1-G 
Tend of blow = 0.9785 x 523°R = 511°R 

From Figure 1-G, actual temperature of the storage tank 
air at the end of a 10 minute run « 510.5*^R. 
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